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Abstract. We present new data structures for approximately counting 
the number of points in an orthogonal range. There is a deterministic 
linear space data structure that supports updates in 0(1) time and ap- 
proximates the number of elements in a 1-D range up to an additive term 
k^^'^ in 0(loglogf/ ■ log log n) time, where k is the number of elements 
in the answer, U is the size of the universe and c is an arbitrary fixed 
constant. We can estimate the number of points in a two-dimensional or- 
thogonal range up to an additive term k'' in 0(log log U + (1/p) log log n) 
time for any p > 0. We can estimate the number of points in a three- 
dimensional orthogonal range up to an additive term k'' in O (log log U + 
(log log 7i)^ + (3" ) log log n) time for « = log ^/ log | -f 2. 

1 Introduction 

Range reporting and range counting are two variants of the range searching 
problem. In the range counting problem, the data structure returns the number 
of points in an arbitrary query range. In the range reporting problem the data 
structure reports all points in the query range. Both variants were studied ex- 
tensively and in many cases we know the matching upper and lower bounds for 
those problems for dimension d < 4. Answering an orthogonal range counting 
query takes more time than answering the orthogonal range reporting query in 
the same dimension. This gap cannot be closed because of the lower bounds for 
the range counting queries: while range reporting queries can be answered in 
constant time in one dimension and in almost-constant time in two and three 
dimensions (if the universe size is not too bigj^, range counting queries take 
super-constant time in one dimension and poly-logarithmic time in two and 
three dimensions. 

Approximate range counting queries help us bridge the gap between range 
reporting and counting: instead of exactly counting the number of points (el- 
ements) in the query range, the data structure provides a good estimation. 
There are data structures that approximate the number of points in a one- 
dimensional interval [4|19j or in a halfspace [7], [15], [2], [8] up to a constant 



^ For simplicity, we consider only emptiness queries. In other words, we ignore the 
time needed to output the points in the answer: if range reporting data structure 
supports queries in 0(/(n) -I- k) time, we simply say that the query time is 0{f{n)). 



factor: given a query Q, the data structure returns the number k' such that 
(1 — e)k < fc' < (1 + e)k, where k is the exact number of points in the an- 
swer and e is an arbitrarily smaU positive constant. In this paper we consider 
the foUowing new variant of approximate range counting: If k is the number 
of points in the answer, the answer to a query Q is an integer fc' such that 
fc — efc" < fc' < fc + efc" for some constant a < 1. Thus we obtain better esti- 
mation for the number of points in the answer for large (superconstant) values 
of fc. On the other hand, if the range Q is empty, then fc' = 0. We present data 
structures that approximate the number of points in a d-dimensional orthogonal 
range for d = 2, 3. We also describe a dynamic one-dimensional data structure. 
Dynamic 1-D Data Structure. A static data structure that answers 1-D re- 
porting queries in 0(1) time is described in [5]. In f5| the authors also describe 
a static data structure that approximates the number of points in a 1-D range 
up to an arbitrary constant factor in constant time. Patra§cu and Demaine |24j 
show that any dynamic data structure with polylogarithmic update time needs 
/2(logn/loglog?7) time to answer an exact range counting query; henceforth 
U denotes the size of the universe. The dynamic randomized data structure of 
Mortensen fl9' supports approximate range counting queries in 0(1) time and 
updates in 0(log'^ U) time; see P2] for other trade-offs between query and up- 
date times. In this paper we present a new result on approximate range counting 
in 1-D: 

— There is a deterministic data structure that can answer one-dimensional 
approximate range counting queries using the best known data structure 
for predecessor queries, i.e. dynamic data structure supports range report- 
ing queries in 0{dpred{n,U)) time, where dpred{n,U) is the time to an- 
swer a predecessor query in the dynamic setting; currently dpred{n, U) = 
0(min(loglogJ7 ■ log log n, ^logn/ log log n)) [6]. We show that we can ap- 
proximate the number of points in the query range up to an additive factor 
fc^/'^, where fc is the number of points in the answer and c is an arbitrary 
constant, in 0{dpred{n, U)) time. We thus significantly improve the preci- 
sion of the estimation; the query time is still much less than the lower bound 
for the exact counting queries in the dynamic scenario. 

Using the standard techniques, we can extend the results for one-dimensional ap- 
proximate range counting to an arbitrary constant dimension d. There is a data 
structure that approximates the number of points in a d-dimensional range up 
to an additive term k'^ for any c > in O(loglogn(logn/loglogn)'^~-'^) time and 
supports updates in 0{log'^^^^'^ n) time. For comparison, the fastest known dy- 
namic data structure [15] supports emptiness queries in O((logn/loglogn)''~^) 
time. Dynamic data structures are described in section [21 

Approximate Range Counting in 2-D and 3-D. We match or almost match 
the best upper bounds for 2-D and 3-D emptiness queries. Best data structures 
for exact range counting in 2-D and 3-D support queries in 0(logn/loglogn) 
and 0((logn/loglogn)^) time respectively [l4] . 
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— If all point coordinates do not exceed n, we can approximate the number of 
points in a two-dimensional query rectangle up to an additive term for 
an arbitrary parameter p, < p < 1, in 0{{1/ p) log log n) time. 

— If all point coordinates do not exceed n, we can approximate the number 
of points in three-dimensional query rectangle up to an additive term in 
0((loglogn)^ -I- (3") log log n) time for an arbitrary parameter p, < p < 1, 
and V = log i/ log | + 2. 

The parameter p is not fixed in advance, i.e. the same data structures can be used 
for answering queries with arbitrary precision. If point coordinates are arbitrary 
integers, then the query time of the above data structures increases by an additive 
term 0(min(log log C/, -^/log n / log log n) ). Data structure for range counting in 
2-D and 3-D are described in section [31 In section [XT] we describe space-efficient 
variants of two- and three-dimensional data structures that estimate the number 
of points in a range up to an additive error fc^ for some fixed constant c. 

Our results for approximate range counting queries are valid in the word 
RAM model. Throughout this paper e denotes an arbitrarily small constant. 

2 Dynamic Approximate Range Counting 

We show that in the dynamic scenario answering one-dimensional counting 
queries with an additive error fc^/"^ can be performed as efficiently as answer- 
ing predecessor queries. The best known deterministic data structure supports 
one-dimensional emptiness queries in 0(dpred(n, U)) time, where dpred{n, U) = 
min ( y/ log ri/ log log n , log log U ■ log log n) is the time needed to answer a prede- 
cessor query in dynamic scenario [S], [Hj. 

Theorem 1. For any fixed constant c > 1, there exists a linear space data struc- 
ture that supports approximate range counting gueries with additive error k^^'^ 
in 0{dpred{n,U)) time, deletions mO(loglogn) amortized time, and insertions 
in 0{dpred{n,U)) amortized time. 

Proof: First we observe that if the query interval contains less than (log log n)'^ 
points for an arbitrary constant c. A; = |P n [a,b]\ < (log log n)'^, then we can 
use a simple modification of the standard binary tree solution: the set P is 
divided into groups of (loglogn)'^ consecutive elements, i.e., \Gi\ = (loglogn)'^ 
and every element in Gi is smaller than any element in G^+i. Using a dynamic 
data structure for predecessor queries we can find in 0{{dpred{n, U)) time the 
successor a' of a in P and the predecessor h' of h in P. If a and h belong to the 
same group Gi, then we can count elements in [a, b] in O(logloglogn) time using 
the standard binary range tree solution. If a' and h' belong to two consecutive 
groups Gi and Gi+i, then we count the number of elements e G Gi, e > a, and 
the number of elements e' G G^+i, e' < b. If a' belongs to a group Gi and b' 
belongs to a group Gj so that j > i + 1, then [a, b] contains more than (log log n)'^ 
elements. We also assume w.l.o.g. that c > 2. 

We maintain the exponential tree [5], [6] for the set P. The root node 
has 0{n^^'^) children, so that each child node contains between n('^~^^/^/2 and 
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2,^(c i)/c points from P. In a general case, if a node v contains points of P, 
then node v has 0{rhj^) children, so that each child contains between n^^ ^^^^ /2 

(c— l)/c 

and 2ny points from P. The exponential tree can be maintained as de- 

scribed in fS], so that insertions and deletions are supported in O (log log n) 
time. Additionally in every node v we store the approximate number of elements 
in any consecutive sequence of children of v, denoted by Ct,(i,j): for any i < j, 

I" 12 < Cv{i,j) < Hy^ + + . . . + Uy. + nj'' 12. 

When r?J'^ 12 elements are inserted into a node v or deleted from w, we set 
Cvihj) = TT-y. + n.„;^j + . . . + iiy. for all i < j. Recomputing Cy{i,j) for a node v 

takes ©(n^^*^) time. Since insertion or deletion results in incrementing or decre- 
menting the value of Uy in O(loglogn) nodes v, recomputing Cy{i,j) incurs 
an amortized cost O(loglogri). Thus amortized cost of a delete operation is 
O(loglogn). When we insert a new point, we also have to find its position in the 
exponential tree; therefore an insertion takes 0{dpred{n, U)) time. 

We store 0(nt^'^) auxiliary values in each node v; hence, we can show that 
the space usage is 0{n) in exactly the same way as in |5l6j . 

Given an interval [a,b], we find b' = pred(6, P) and a' = succ(a,P) and 
identify the leaves of the exponential tree in which they are stored. The lowest 
common ancestor q of those leaves can be found in O (log log n) time because 
the height of the tree is O (log log n). If a' and 6' are stored in the i-th and 
the j-th children of q and i + 1 < j, then all elements stored in qi+i, . . . , qj-i 
belong to [a,b] and we initialize a variable count to Cv{i + 1, j — I) ■ Otherwise 
count is set to 0. Then, we traverse the path from q to a' and in every visited 
node V we increment count by Cy(iv + l,ry), such that a' is in the iy-th child 
of V, and r„ is the total number of w's children. Finally, we traverse the path 
from q to b' and in every visited node v we increment count by Cy{l,iy — 1), 
such that b' is in the iy-th child of v, Suppose that the variable count was 
incremented by s„ > when a node v was visited. Let ky be the exact number 
of elements in all children of v whose ranges are entirely contained in v. Then, 

ky TZy ^ Sy ^ ky ~\~ny . SlUCC ky ^ , ky ky ^ Sy ^ ky~\~ky 

Clearly, the total number of points equals to the sum of ky for all visited nodes v. 
The search procedure visits less than Ch log log n nodes for a constant Ch- Hence, 
k — loglogn < count < k + fc3/('=~i) loglogn for k = \P H [a, b]\. Since 

loglogn < fci/^'^-i), k - < count < k + k^/^''-'^\ We obtain the result 

of the Theorem by replacing c with c' — max(5c, 5) in the above proof. □ 

Our dynamic data structure can be extended to d dimensions using the standard 
range tree [TO]. 

Theorem 2. For any fixed constant c > 1, there exists a data structure that sup- 
ports d- dimensional approximate range counting queries with additive error k^^'^ 
in 0(log log n(logn/ loglogn)''"^) time and updates in OiXog'^^^^'^ n) amortized 
time. 
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Proof: This result can be obtained by combining the standard range tree tech- 
nique (node degree in a range tree is 0(log'^ n) for an appropriate constant 
e' = e/{d — 1)) with the data structure for one-dimensional approximate range 
counting of Theorem[TJ Details will be given in the full version of this paper. □ 

3 Approximate Range Counting in 2-D and 3-D 

A point p dominates a point q if each coordinate of p is greater than or equal 
to the corresponding coordinate of q. The goal of the (approximate) dominance 
counting query is to (approximately) count the number of points in P that 
dominate q. The dominance query is equivalent to the orthogonal range query 
with a restriction that query range Q is a product of half-open intervals. We start 
this section with a description of the data structure that estimates the number 
of points in the answer to a 2-D dominance query up to a constant factor. We 
can obtain a data structure for general orthogonal range counting queries using 
a standard technique. Then, we show that queries can be answered with higher 
precision without increasing the query time. Finally, we describe a data structure 
for approximate range counting in 3-D. For simplicity, we only consider the case 
when all point coordinates are bounded by n. We can obtain the results for the 
case of arbitrarily large point coordinates by a standard reduction to rank space 
technique 13J: the space usage remains linear and the query time increases by 
pred{n, U) - the time needed to answer a static predecessor query. 

Theorem 3. There exists a linear space data structure that answers approx- 
imate two-dimensional dominance range counting queries on n x n grid in 
O(loglogn) time. 

A ^-approximate boundary, introduced by Vengroff and Vitter ^26j is a polyline 
M consisting of 0{n/t) axis-parallel segments that partitions the spacci, so 
that every point M is dominated by at most 2t and at least t points of P. This 
notion can be straightforwardly extended to a ta-boundary tVJq,: M.a partitions 
the space into two parts, and every point A^^ is dominated by at most a ■ t 
and at least t points of P. We can construct a i^-boundary with the same 
algorithm as in [55]. Let p be a point with coordinates (0,0). We move p in 
the positive x direction until p is dominated by at most at points. Then, we 
repeat the following steps until the x-coordinate of p equals to 0: a) move p in 
-\-y direction as long as p is dominated by more than t points of P b) move p 
in the —x direction until p is dominated by at points of P. The path traced 
by p is a ia-boundary; see Fig. [1] for an example. Inward corners are formed 
when we move p in -\-y direction, i.e. inward corners mark the beginning of 
step a) resp. the end of step b). Inward corners of J\A have a property that 
no point of Ai is strictly dominated by an inward corner and for every point 
m G that is not an inward corner, there is an inward corner dominated by 
m. There are 0{n/t) inward corners in a ^Q-approximate boundary because for 

^ In this section we assume that all points have positive coordinates 
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Fig. 1. Example of f-approxiinate boundaries in 2-D. For simplicity, the points 
of the set P are not shown. 

every inward corner c = {c^, Cy) there are (a — l)t points that dominate c and 
do not dominate inward corners whose a;-coordinates are larger than c^- Our 
data structure consists of log^ n tct-approximatc boundaries M.1.M21 ■ ■ ■ 
such that M-i is an a*-approximate boundary of P, i.e. every point on M.i is 
dominated by at least a' and at most a'+^ points of P. If a point p € Mi is 
dominated by a query point then q is dominated by at most a^~^^ points of P. 
If q dominates a point on M.i, then it also dominates an inward corner of A^j. 
Hence, we can estimate the number of points that dominate up to a constant 
a by finding the minimal index j such that q dominates an inward corner of M.j . 
Since q is dominated by a point of M.j-i, q is dominated by A: > a^~^ points of 
P. On the other hand, k < a^^^ because a point of M.j is dominated by q. 

Wc can store inward corners of all boundaries M.i in a linear space data 
structure so that for any point q the minimal index j, such that some point 
on M.j is dominated by q, can be found in O(loglogn) time. We denote by 
prcd^(a, S") the point p = {px,Py) G S, such that px = prcd(a, S'^;) where Sx 
is the set of x-coordinates of all points in S. For simplicity, we sometimes do 
not distinguish between a boundary Mi and the set of its inward corners. Let 
q = {qx,<ly)- Let Cj = {cx,Cy) be the inward corner on a boimdary Aii whose 
x-coordinate precedes qx, Ci = pred^(ga:, Alj). For any other inward corner 
= (c^,Cy) on A^j, c'y > Cy if and only if < because the y-coordinates of 
inward corners decrease monotonously as their x-coordinates increase. Hence, q 
dominates a point on A4i if and only if qy > Cy. Thus given a query point q, it 
suffices to identify the minimal index j, such that the y-coordinate of the inward 
corner Cj G A4j that precedes qx is smaller than or equal to qy. The .T-axis is 
subdivided into intervals of size logn. For each interval Is the list Lg contains 
indexes of boundaries Mi such that the x-coordinate of at least one inward corner 
of Mi belongs to Ig. For a query point q with qx G Ig and for every j G Lg, 
we can find the inward corner preceding qx with respect to its x-coordinate, 
pTedx{qx,Mj), in 0(1) time because x-coordinates of all relevant inward corners 
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belong to an interval of size logn. Hence, we can find the minimal index js S Ls, 
such that q dominates a point on A4j^ in O(loglogn) time by binary search 
among indexes in Lg. For the left bound of an interval Ig = [as,bs] and for 
all indexes j = 1, . . . , log„ n, the list Ag contains the inward corner Cj, such that 
Cj = pTed^{as,Mj). By binary search in we can find the minimal ja such that 
q dominates the inward corner Cj^ £ As- Clearly j = mm{ja,js) is the minimal 
index of a boundary dominated by q. 

Theorem 4. There exists a O(nlog^n) space data structure that supports two- 
dimensional approximate range counting queries on n x n grid in O(loglogn) 
time. 

The next Lemma will enable us to obtain a better estimation of the number of 
points. 

Lemma 1. There exists a 0{n\ogn) space data structure that supports two- 
dimensional approximate range counting queries on n x n grid with an additive 
error in 0((l/p) log logn) time for any p, < p < 1. 

Proof: We divide the grid into x-slabs Xi = [xi-i,Xi] x and y-slabs Yj = 
[1, n] X [Vj-i, Vj], so that each slab contains n^^^ points. For every point (.7:,;. jjj), 
< < n^/^ we store the number of points in P that dominate it. There is 
also a recursively defined data structure for each slab. The total space usage is 
s{n) = 0{n) + 2n^/^s(n^/^) and s{n) = O(nlogn). 

We can easily obtain an approximation with additive error 2n^/^ using the 
first level data structure: for a query q = {qx,qy) we identify the indexes i and j, 
such that Xi-i < (/x < Xi and yj-i < qy < yj, i.e. we identify the .x-slab Xi and 
the y-slab Yj that contain q. Indexes i and j can be found in 0(log log n) time. Let 
c{x, y) be the number of points that dominate a point p = {x, y); let c{x, y, Xi) 
{c{x, y, Yj)) be the number of points in the slab Xi (Yj) that dominate p = (x, y). 
Then c{qx, qy) = c{xi, yj) + c{xi, qy, Yj) + c{qx,qy, Xi). Since c{xi, qy, Yj) < n^l'^ 
and c{qx,qy,Xi) < n^/^, the value of c{xi,yj) is an approximation of c{q,^,qy) 
with an additive error 2ni/2. Using recursive data structures for slabs Xi and 
Yj we can estimate c{qx,qy, Xi) and c{xi,qy,Yj) with an additive error 2n^/'^ 
and estimate c{qx,qy) with an additive error 4n^/^. If the recursion depth is v 
(i.e. if we apply recursion v times), then the total number of recursive calls is 
0(2'') and we obtain in 0((2'') log logn) time an approximation with additive 
error 2" ■ n^/^ for any positive integer v. 

We set recursion depth v = [log(l/p)] + 2. Then, v + (1/2") logn < 
(p/4)logn + log(l/p) = (p/4+ i2^M)iogn < plogn. Hence, n^ > 2''n^/^\ 
Therefore, if recursion depth is set to v, then our data structure provides an 
answer with additive error n''. □ 

Theorem 5. There exists a 0(n log^ n) space data structure that supports two- 
dimensional dominance counting queries on n x n grid with an additive error k'' 
for an arbitrary parameter p, < p < 1, in 0{{l/p) log logn) time. 
There exists a 0(n log^n) space data structure that supports two-dimensional 
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range counting queries on n x n grid with an additive error for an arbitrary 
parameter p, Q < p < 1, in 0((l//3) loglogn) time. 

Proof: As in Theorem [3] we construct ^-boundaries Ali, . . . , A^iogn, such that 
Mi is a 2'-approximate boundary, i.e. each point on M.^ is dominated by at least 
2' and at most 2^' points of P. For each inward corner Cij of every Mj, we store 
a data structure Dij that contains aU points that dominate Ci^j and supports 
approximate counting queries as described in Lemma [TJ For a fixed j , there are 
0(|j) data structures Di j, and each Di j contains 0(2^) points. Hence, all data 
structures Dij use O(nlog^n) space. 

As described in Theorem|4l we can find in 0(log log n) time the minimal index 
j, such that Aij is dominated by the query point q and an inward corner a^j G 
Mj dominated by q. Then, we use the data structure Di^j to obtain a better 
approximation. Since Dij contains 0{k) points, by Lemma[l]Z?i.j estimates the 
number of points that dominate q with an additive error in 0((l/p) log log n) 
time. We can extend the result for dominance counting to the general three- 
dimensional counting using the standard technique from range reporting |12I25] ; 
see also the proof of Theorem □ 

Lemma 2. There exists a 0{n\o^ n) space data structure that supports three- 
dimensional approximate range counting queries onnxnxn grid with an additive 
error in 0{3^' log log ri) time for any p, < p < 1, and for v = log log § + 2. 

Proof: We divide the grid into x-, y-, and z-slabs, Xi — [2:^-1, Xi] x [1, n] x [1, n], 
Yj — [l,n] X [yj-i,yj\ x [\,n], Zd — [^,n] x [l,n] x [zd-i,Zd], so that each slab 
contains n^/^ points. For each point {xi,yj, Zd) we store the number of points in 
P that dominate it. There is also a recursively defined data structure for each 
slab. The total space usage is s(n) = 0{n)-\-'3n^^^s{n'^/^) and s(n) — O(nlog^n). 

For a query q = {q^, qy,qz) we identify the x-, y-, and z-slabs Xi, Yj, and Zd 
that contain q. By the same argument as in Lemma[Tl the number of points that 
dominate {xi,yj,Zd) differs from the number of points that dominate q by at 
most 3n^/'^. We can estimate the number of points that dominate q and belong 
to one of the slabs Xi, Yj, and Zd using recursively defined data structures. If 
the recursion depth is v, then we obtain in 0(3" log log n) time an approximation 
with additive error 3" •n'^^/'^^ for any positive integer v. The result of the Lemma 
follows if we set v — log i/ log | + 2. □ 

Theorem 6. There exists a 0(nlog'*n) space data structure that supports ap- 
proximate dominance range counting queries on n x n x n grid with an additive 
error k^in 0((loglogn)^ -|- 3^1oglogn) time for any p, < p < \, and for 
V = log i/ log I -I- 2. 

There exists a O(nlog^n) space data structure that supports approximate range 
counting queries onnxnxn grid with an additive error k^ in 0((loglogn)^ -|- 
3" log log n) time for any p, < p < 1, and for v — log ^/ log § + 2. 
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Proof: Instead of counting points that dominate q we count points dominated 
by q. Both types of queries are equivalent. Hence, the data structure of Lemma[2] 
can be used to approximately count points dominated by q. 

A downward corner of a point p consists of all points dominated by p. We 
define an approximate t-level as a set of downward corners C, such that (1) 
any point p that dominates at most t points of P is contained in some r £ C 
(2) any downward corner r £ C contains at most a ■ t points of P. Afshani [T] 
showed that for an arbitrary constant a there exists an approximate t-level of 
size 0{j). We can assume that no r e £ dominates r' £ C va an approximate 
t-level C: if r dominates r', then the downward corner r' can be removed from C. 
Identifying an inward corner r £ C that dominates a query point q (or answering 
that no r G £ dominates q) is equivalent to answering a point location query in 
a rectangular planar subdivision 26 21j and takes 0((loglogn)^) time. 

Our data structure consists of approximate levels A^i, • ■ ■ j -^logru such 
that M-i is a 2*-approximate level and the constant a is chosen to be 2. For 
every downward corner r^ j G Alj, we store all points dominated by r^j in a 
data structure Dij] Di j contains 0(2^) points and supports counting queries 
with additive error 0(2^^) by Lemma[2l All data structures Dij use 0(n log'* n) 
space. 

We can find a minimal j, such that A4j dominates q in 0((loglogn)^) time 
by binary search. Let r^j be the downward corner that dominates q. We can 
use the data structure Di j to estimate the number of points that are dominated 
by q with an additive error k^; by Lemma [2] this takes 0(3"loglogn) time for 
^; = logi/log| + 2. 

We can extend the result for dominance counting to the general three- 
dimensional counting using the standard technique |12l25j : see also the proof 
of Theorem H □ 

3.1 Space-Efficient Approximate Range Counting in 2-D and 3-D 

If we are interested in counting with an additive error fc"^ for some predefined 
constant c > 0, then the space usage can be significantly reduced. The two- 
dimensional data structure uses O(nlog^n) space (0(n) space for dominance 
counting), and the three-dimensional data structure uses 0{n log^ n) space (0(n) 
space for dominance counting). The main idea of our improvement is that in 
the construction of Lemma [1] (resp. Lemma [2]) each slab contains n^^^'^'^ points 
^j^2/3-i-e points) for some e > and there is a constant number of recursion levels. 

Lemma 3. For any fixed constant c < 1, there exists a 0(rt*^^) space data 
structure that supports two-dimensional approximate range counting queries on 
n X n grid with an additive error n'^ in O (log log n) time. 

Proof: We divide the grid into x-slabs Xi = [xi^i,Xi] x [l,n] and y-slabs Yj = 
[l,n] X [yj_i,?/j], so that each slab contains points. As in Lemma[l] we 

store for each point {xi, yj), < i,j, < n*/^~'^, the number of points in P that 
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dominate it. Note that there are 0{n^~^^) points {xi,yj) for <i,j < rt^/^^^ . 
If an x-slab or a y-slab contains more than n-^ points for a constant / = c/4, we 
store a recursively defined data structure for that slab. The number of recursion 
levels is g = r iQg(2/^(i4^2e)) 1 • ^i^*^'^ each point is stored in one recursively defined 
data structure for an x-slab and in one recursively defined data structure for 
a y-slab, the total number of points in all recursively defined data structures 
increases by factor 2 with each recursion level. Thus the total space usage is 
ELi2^-0(ni-) = 0(ni-). 

Given a query q = (qx^Qy), we identify the x-slab Xi and the y-slab Yj that 
contain q. Let c{x, y) be the number of points that dominate a point p = (x, y); let 
c(a;, y, Xi) {c{x, y, Yj)) be the number of points in the slab Xi {Yj) that dominate 
p = {x,y). As in the proof of Lemma [1] c{qx,qy) = c(xi,yj) + c(xi,qy,Yj) + 
c{qa:,qy, Xi)^ where Xi and Yj are the z-slab and the y-slab that contain q. 
If slabs Xi and Yj, contain more than points, we estimate c{xi,qy,Yj) and 
c{qx, qy, Xi) using data structures for slabs Yj and Xi. Otherwise we use c{xi, yj) 
as an estimation for c{qx, qy). By the same argument as in the proof of Lemma[TJ 
we obtain an approximation with additive error 2^ -n^ . Since g < 21og(l//) and 
/ = c/4, g + /logn < 21og(l//) + (c/4) logn < clogn. Hence, 2^ ■ < n'^ and 
we estimate the number of points in a range with an additive error that is less 
than n'^. □ 

Using Lemma [21 we can prove the following Theorem. 

Theorem 7. For any fixed constant c < 1, there exists a 0(n) space data struc- 
ture that supports two-dimensional dominance counting queries on n x n grid 
with an additive error k'^ in O(loglogn) time. 

For any fixed constant c < 1, there exists a 0{n\o^ n) space data structure that 
supports two-dimensional range counting queries on n x n grid with an additive 
error k'^ mO(loglogn) time. 

Proof: We construct a sequence of t-approximate boundaries A4i in the same 
way as in Theorem [5] and store all points that dominate an inward corner aj in 
data structure Dij. The only difference is that Dij is implemented as described 
in LemmaO For a fixed j, there are 0{~) data structures Dij, and each Dij 
needs 0{2'^^~'^^'^) space. Hence, all data structures Dij use 0{^j jft) = 0{n) 
space. 

Dominance queries are processed in exactly the same way as in Theorem [5l 
We can extend the result for dominance counting to the general two-dimensional 
counting using the standard technique from range reporting [12|25j ; see also the 
proof of Theorem m □ 

Lemma 4. For any fixed constant c < 1, there exists a 0{n^^^) space data 
structure that supports three-dimensional approximate range counting queries 
on n X n X n grid with an additive error n'^ in O(loglogn) time. 

Proof Sketch: Like in Lemma [2l we divide the grid into x-, y-, and z-slabs, 
X, = [x^-l,x^] X [l,n] X [l,n], Yj = [l,n] x [yj_i,yj] x [l,n], Za = [l,n] x [l,n] x 
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[zd-i,Zd], but each slab contains n^/^+^ points. For each point {xi,yj,Zd) we 
store the number of points in P that dominate it. If the number of points in a 
slab is greater than n-^ for / = c/16, then we store a recursively defined data 
structure for each slab. 

We can estimate the space usage and analyze the query algorithm in the 
same way as in Lemma [3l □ 

Theorem 8. For any fixed constant c < 1, there exists a 0{n) space data struc- 
ture that supports approximate dominance range counting queries on n x n x n 
grid with an additive error k'^ in 0((loglogn)'^) time. 

For any fixed constant c < 1, there exists a 0{n log^ n) space data structure that 
supports approximate range counting queries on n x n x n grid with an additive 
error k'^ m 0((loglogn)'^) time. 

Proof Sketch: As in the proof of Theorem [6] our data structure consists of 2*- 
approximate levels Aii for i — 1, . . . ,logn. For every inward corner r^j G M.j, 
we store all points dominated by ri.j in the data structure described in 
LemmaHl Each Dij uses 0(2^^"^^^ ) space. Since a 2^ -approximate level M.j has 
0{~) inward corners, all Mj use Oi^,- — 0{n) space. 

Dominance counting queries are answered in the same way as in Theo- 
rem |6l We can extend the result for dominance counting to the general three- 
dimensional counting by applying the standard technique from range report- 
ing |12|25j that was also used in proofs of Theorems [H [6l [71 □ 
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Appendix A. Proof of Theorem |4] 

We use the well known technique used for range reporting queries |12I25| . The 
set of points P is subdivided into subsets Pi , P2 , • ■ • , -Ps , so that the total number 
of points in Pi U . . . U is 0(nlog^ n), and an arbitrary query rectangle Q can 
be represented as a union of at most four rectangles Qi, . . . , Qs, s < 4, so that 
Q n P = (Qi n Pii) U . . . U [Qs n P^J and each Qi is a product of two half- 
open intervals. We store the date structure for approximate dominance queries 
of Theorem [3] for each set Pi, so that the total space usage is 0(n log^ n). Given 
a query Q, we can decompose Q into Qi, ■ ■ ■ ,Qs and find the corresponding 
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Pi-^, . . . ,Pi^ in O(loglogn) time, see e.g. ^21]. Then, we can estimate the number 
of points in each Pi CiQj, 1 < j < s, and thus estimate the number of points in 

Q n P = (Qi n j'u . . . u (Q, n J 
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